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Chapter 25 

Quantum Phase Transitions 
in Dense QCD 

Quantum chromodynamics (QCD) at finite temperature, T, and quark chemi- 
cal potential, /x, has a rich phase structure: at low T and low /i, the Nambu- 
Goldstone (NG) phase with nearly massless pions is realized by the dynamical 
breaking of chiral symmetry through condensation of quark— anti-quark pairs, 
while, at low T and high /i, a Fermi liquid of deconfined quarks is expected to 
appear as a consequence of asymptotic freedom. Furthermore, in such a cold 
quark matter, condensation of quark— quark pairs leads to the color supercon- 
ductivity (CSC). At high T for arbitrary fj,, all the condensates melt away and 
a quark-gluon plasma (QGP) is realized. The experimental exploration of ther- 
mal phase transition from the NG phase to QGP is being actively pursued in 
ultrarelativistic heavy ion collisions at RHIC (Relativistic Heavy Ion Collider), 
and will be continued in the future at LHC (Large Hadron Collider). The quan- 
tum phase transition from the NG phase to the CSC at low T is also relevant 
to heavy-ion collisions at moderate energies, and is of interest in the interiors 
of neutron stars and possible quark stars. 

In this Chapter, after a brief introduction to the basic properties of QCD, 
the current status of the QCD phase structure and associated quantum phase 
transitions will be summarized with particular emphasis on the symmetry re- 
alization of each phase. Possible connection between the physics of QCD and 
that of ultracold atoms is also discussed. 

25.1 Introduction to QCD 

The color SU(3)c gauge theory of quarks and gluons [1] is now called the quan- 
tum chromodynamics (QCD) and is established as the fundamental theory of 
strong interaction. The Lagrangian density of QCD reads 

£qcd = 9 L ^9 L + q R Hpq R - \g%G^ + q L mq R + g R m<? L , (25.1.1) 
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Figure 25.1: The QCD fine-structure constant a s determined from the r decay, 
the T decay, the deep inelastic scattering, the e + e~ annihilation and the Z- 
boson resonance shape and width [3]. 



where the covariant derivative is defined as = 4- igt™A%, with g being 
the QCD coupling constant, t a the SU(3)c group generator and A£ the gluon 
field belonging to the adjoint representation of SU(3)c- The gluon field-strength 
tensor is defined as <3g" = 9M£ - d v A£ - gf^A^A^ with f a ^ being the 
structure constant of the SU(3)c group. 

The quark field q belongs to the fundamental representation of SU(3)c- The 
right (left) handed quark q n = |(1 + j 5 )q (q L = |(1 — 75)9) is an eigenstate of 
the chirality operator 75 with the eigenvalue +1 (—1). Although quarks have 
six flavors (u,d,c,s,t,b) in the real world [2J, we focus only on three light quarks 
(u,d,s) in this Chapter, so that the quark mass matrix is m = diag(m u , rrid, m s ). 
As is evident from Eq. (|25.1.1[) . only the mass term can mix the left-handed quark 
and the right-handed quark in the QCD Lagrangian. 

The running coupling constant g(n) is defined as an effective coupling strength 
at the energy scale k. Due to the asymptotic freedom of QCD [3], g(n) becomes 
small when k increases as seen explicitly in the two-loop perturbation, 



a s («) 



g 2 («) 

4.7T 



W ln(^/A2 J 



Pl ln(ln(^/A^ c J) 
P 2 o ln(«7A2 CD ) 



(25.1.2) 



f Af/)/(47r) 4 with Nf being 



where f3 = (11 - §A7)/(4tt) 2 and /3 X = (102 

the number of flavors. Here, A QCD (~ 200 MeV = 2 x 10 s eV) is called the 
QCD scale parameter which is determined from the comparison of Eq. (|25.1.2[) 
with the experimental data in high energy processes satisfying k 3> A QOD (sec 

Figia. 

Equation (|25.1.2|) implies that a s (n ~ A QCD ) ~ O(l), so that the QCD per- 
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Figure 25.2: The masses of u, d and s quarks at the scale k = 2 GcV = 2 x 10 9 eV 
determined from various observables and methods [4]. 



turbation theory breaks down. This leads to various non-perturbative phenom- 
ena such as the confinement of quarks and gluons and the dynamical breaking 
of chiral symmetry [51 [5] at low energies, k < A QCD . These are responsible for 
the formation of hadrons (qq mesons and qqq baryons), and also for the ori- 
gin of their masses. On the other hand, at extremely high temperature and/or 
high baryon density where ct s (k S> A qcd ) < 1, the system may be treated as a 
weakly interacting matter of quarks and gluons. Thus, there must be a phase 
transition from the hadronic matter composed of confined quarks and gluons at 
low energies to the deconfined quark-gluon matter at high energies. 

Due to quantum corrections, the quark mass m also becomes k dependent. 
As seen from Fig |25.2l the current determination of the u and d quark masses at 
K = 2 GeV indicates that they are about 50 to 100 smaller than A pCD , while s 
quark mass is comparable to A QCD . Therefore, it is legitimate to treat m u /A pCD 
and rrid/ A QCD as small expansion parameters, while the expansion by m s /A pCD 
does not necessarily work. Systematic expansion in terms of the quark masses 
is called the chiral perturbation theory and has been successfully applied to a 
wide variety of QCD phenomena [7]. 



25.1.1 Symmetries in QCD 

Let us consider the following transformations of the quark fields, 



e ~*e B e - t e A Vl Vc q ^ 



H8 A 



V n V c q R , (25.1.3) 



where Vc (gauge rotation) is a local SU(3)c transformation in the color space, 
while Vl(r) (chiral rotation) is a global SU(3)l(r.) transformation in the flavor 
space. The 9b and 8a are phases associated with a global U(1)b transformation 
(baryon-number rotation) and the global U(1)a transformation (axial rotation), 
respectively. For m u .d, s = (the flavor-SU(3) chiral limit), QCD Lagrangian 
Eq. (|25.1.1[) is invariant under Eq. (|25.1.3[) together with the SU(3)c gauge trans- 
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Figure 25.3: Light hadron spectra obtained from lattice QCD Monte Carlo 
simulations with dynamical u, d, s quarks in the Wilson fcrmion formalism. The 
spatial lattice volume V and the lattice spacing a are (2.9 fm) 3 and 0.09 fm, 
respectively. Horizontal bars denote the experimental values |10) . 



formation of the gluons, so that the full continuous symmetries of QCD become 



= [SU(3) c ] local ® [SU(3) L ® SU(3) R ] global ® [U(l) B ] global . (25.1.4) 

Although the U(1)a looks like a symmetry of Eq. (|25.1.1j) . it is explicitly broken 
by quantum effect known as the axial anomaly [8J which reduces U(1)a down to 
its discrete subgroup Z(27V/)a = Z(6)a- The masses of light quarks m U! d, s 
act as small external fields to break the global chiral symmetry [SU(3) L (g> 

SU (3)rW 

In the past few years, there arises a remarkable progress in lattice gauge 
theory |5] particularly in calculating the hadron spectra on the basis of the 
QCD Monte Carlo simulations with light dynamical u, d, s quarks. This has 
been achieved partly due to the growth of the supercomputer speed and partly 
due to the new algorithms: Simulations with quark masses very close to the 
physical values are now possible in the Wilson fermion formalism [10[lllj. Shown 
in Fig l25.3l is an example of such calculations for meson and baryon masses 
extrapolated to the physical quark masses using the simulation data taken in 
the interval, ^(m u + md) = 3.5 MeV — 67 MeV at k — 2 GeV. The experimental 
values are reproduced in 3% accuracy at present. 
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25.1.2 Dynamical breaking of chiral symmetry 

Although QCD Lagrangian in the flavor-SU(3) chiral limit has the symmetry Q 
in Eq. (|25. 1 .4[) , the ground state of the system breaks some of the symmetries dy- 
namically. Consider the QCD vacuum |0) at zero temperature and zero baryon 
density. Taking into account the fact that the QCD does not allow dynamical 
breaking of parity and vector symmetries in the vacuum |12j . the following is 
one of the possible symmetry breaking patterns, 

g -» SU(3) C ® SU(3) L+R ® U(1) B , (25.1.5) 

where simultaneous transformation of the left-handed and right-handed quarks 
(the vector rotation, Vl = Vr) as indicated by SU(3) L+R remains as a symmetry 
of the vacuum. An order parameter to characterize this "dynamical breaking of 
chiral symmetry" would be the "chiral condensate", (qq)o = (0|gg|0), which is 
not invariant under the opposite rotation of the left-handed and right-handed 
quarks, Vl = Vp\ Recent lattice QCD simulation using overlap Dirac fermion 
with dynamical u, d, s quarks indeed shows that Eq. (|25.1.5|) is realized with the 
chiral condensate [13] 

i(uu + dd) ^ -(242(04) MeV) 3 at k = 2 GeV. (25.1.6) 

The non-vanishing chiral condensate, (qq)o — (q^q^ + <7 R 9 L )o ^ 0' implies that 
the quark— anti-quark pairs are Bose-Einstein condensed. It also implies that a 
non-perturbative mixing between the left-handed and the right-handed quarks 
takes place in the QCD vacuum: In other words, an effective quark mass called 
the chiral gap is dynamically generated. Indeed, there are phenomenological 
evidences that the u, d quarks and s quark have effective masses M u d ~ 350 
MeV and M s ~ 550 MeV inside hadrons [5J. 

The Nambu-Goldstone (NG) bosons associated with the dynamical break- 
ing of the flavor-SU(3) chiral symmetry are nothing but the pions, kaons and 
the 77-meson. Moreover, one can derive a spectral sum rule, called the Gell- 
Mann— Oakes— Renner relation [14], which relates the pion mass to the chiral 
condensate as f%m*± — — rh(uu + dd}o + 0(rh 2 ) . Here m = (m u + md)/2, f n ( = 
92.4 MeV) is the pion decay constant, and m^i (~ 140 MeV) are the charged 
pion masses. Similar relation holds also for the neutral pion ir°. In the limit 
w U)( j — * 0, the pion mass vanishes as it should be from the Nambu-Goldstone 
theorem. 

25.2 QCD matter at high temperature 

As the temperature T of the system increases, the condensed qq pairs in the 
QCD vacuum are melted away by thermal fluctuations. This is analogous to the 
phase transition in metallic superconductors with the electron pairing (e-|-ej,) as 
an order parameter. 
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Figure 25.4: Normalized chiral condensate lZq q = [(uu) — (m/m s )(ss)]/[(uu)o — 
(m/m s )(ss)o] as a function of T for two different lattice spacings, a = 0.24 fm 
and a = 0.17 fm, calculated by the lattice QCD simulations with dynamical u, 
d, s quarks in the staggered fcrmion formalism. Vertical band in the middle 
indicates the pseudo-critical temperature T pc |15j . 



To see that (qq) vanishes at extreme high T, let us consider the QCD parti- 



tion function at zero baryon density, Zqcb = Tr 



-Hqco/T 



3 P(T)V/T which 



leads to (qq) = — ■ If we have a situation where T A QCD , the system is 
approximated by the Stefan-Boltzmann gas of free quarks and gluons because of 
the asymptotic freedom. Since the quark-gluon vertex does not change chirality 
(the transition between q L and q R is not allowed in perturbation theory) and 
thus the expectation value of qq = q^q R + <7 R 9 L vanishes in any finite order of 
the perturbation as long as m q = 0. 

As shown in Fig ]25.4i lattice QCD Monte Carlo simulations at finite T with 
zero baryon density indeed indicate a sudden drop of the chiral condensate 
around the pseudo-critical temperature determined from the susceptibility peak 

nsum, 

T pc ~ (150 - 200) MeV = (1.7 - 2.3) x 10 12 K. (25.2.7) 

Note that the phase transition from the low T phase (Nambu-Goldstone phase) 
to the high T phase (quark-gluon plasma) is first order in the flavor-SU(3) 
chiral limit (m U! d, s = 0) and is second order in the flavor-SU(2) chiral limit 
(w u ,d = 0, m s = oo ) [18j . On the other hand, in the presence of the finite 
quark masses m Ui d, s acting as external fields, the transition is crossover as seen 
in Fig l25.4l For more details on the QCD thermodynamics, see the review [19]. 

The high temperature quark-gluon plasma (QGP) is believed to be present 
in the early universe with its age younger than 10 -5 sec. Attempts to create 
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Figure 25.5: Possible internal structures and compositions of four different types 
of compact stars [2U] • 



such extremely hot system by the relativistic heavy-ion collisions in the labo- 
ratory have been started from 2000 at RHIC (Relativistic Heavy Ion Collider) 
in Brookhaven national laboratory and will be pursued at LHC (Large Hadron 
Collider) in CERN. RHIC has already produced a plenty of data showing not 
only the evidence of QGP but also strongly interacting characters of QGP [2D] . 
(See the recent reviews |21) on this rapidly developing subject.) In the following, 
we will focus more on QCD matter with finite baryon density at low T. Such 
a system may undergo successive quantum phase transitions which are relevant 
to the physics of neutron stars and of possible quark stars. 



25.3 QCD matter at high baryon density 

Soon after the discovery of the asymptotic freedom of QCD, a possible transition 
from hadronic matter to quark matter in the core of the neutron stars has been 
pointed out [22 . A strange quark star entirely made of deconfmed u, d, s quarks, 
yet undiscovered, was also proposed |23j as a modern version of the early idea 
of the quark star [21]. Although there has been no observational evidence of 
quark stars yet, nature may be strange enough to accommodate such compact 
object in our universe and we should prepare for the future discovery. 

Shown in Fig l25.5l is a schematic view of various forms of compact stars, from 
the neutron star to the quark star. Typical radius of the neutron star is about 
10 km, while its mass is comparable to the solar mass M (~ 2 x 10 30 kg). 
Since neutrons cannot be bound by the strong interaction only, the presence 
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of gravitational force is essential to hold the neutron star. This implies that 
the radius increases as the mass decreases. The masses of the neutron stars in 
the binary systems are centered around 1.35M0. Observed upper limit of the 
surface temperature of neutron stars is less than 10 9 K after 1 year from their 
births. In the early stage, the cooling occurs through the neutrino emissions, 
while in the later stages it is dominated by surface photon emissions. 

For pulsars (rotating neutron stars), the measured rotational frequency is 
ranged from milli second to several seconds. The surface magnetic field is typ- 
ically 10 12 gauss for ordinary pulsars with rotational period P ~ Is and 
dP/dt ~ 10 -15 . There are also stars with much larger (smaller) dP/dt and 
larger (smaller) magnetic field ~ 10 15 (10 9 ) gauss. Sudden spin up of the rota- 
tion associated with a subsequent relaxation to the normal rotation has been 
observed and is called the glitch. This phenomena should be related to the in- 
ternal structure of neutron stars, in particular the superfluidity of the neutron 
liquid. 

The outer crust of the neutron stars is a solid composed of heavy nuclei 
forming a Coulomb lattice in the sea of degenerate electrons. As the pressure 
and the density increase toward the inner region, electrons tend to be captured 
by nuclei and at the same time neutrons drip out from the nuclei, so that 
the system is composed of neutron-rich heavy nuclei in the Fermi sea of the 
neutrons and electrons. Eventually, the nuclei dissociate into neutron liquid 
and the system becomes a degenerate Fermi system composed of superfluid 
neutrons together with a small fraction of superconducting protons and normal 
electrons. 

When the baryon number density (p) of the core of the neutron stars exceeds 
a few times of the central density of heavy atomic nuclei (po = 0.16 fm~ 3 = 
0.16 x 10 39 cm -3 ), one may expect exotic components such as the hyperons 
(baryons with s quarks), Bose-Einstein condensates of pions and kaons, and 
the deconfined quark matter, which can contribute to the acceleration of the 
neutron-star cooling. For more details of the physics of high density matter and 
compact stars, see e.g. [25| [26] . 

25.3.1 Neutron-star matter and hyperonic matter 

Although the neutron star is mainly composed of the degenerate neutrons, 
other species are also present as a result of the chemical equilibrium condi- 
tions. Indeed, the matter made of only neutrons is unstable against the /3-decay, 
n — > p + e~ + P c . After the decay, the electron- neutrino leaves the star with- 
out much interactions if the neutron star is cold enough. On the other hand, 
the protons and the electrons remain in the star and form degenerate Fermi 
liquid together with the neutrons. The equilibrium configuration of n, p and 
e~, which we call the standard neutron-star matter (see Fig |25.6f a)), is deter- 
mined by the three conditions: chemical equilibrium, charge neutrality and the 
baryon- number conservation: p n = fi p + /x e , p p = p e , p = p n + p p , where pi 
denotes the number density of z-species. 

If we assume the non-interacting degenerate fermions for simplicity, it is 
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Figure 25.6: Compositions of matters under chemical equilibrium and charge 
neutrality conditions in the Fermi gas model. Shaded areas show the occupied 
states [20 . (a) The neutron-star matter with n, p, and e~, (b) the hyperon 
matter with n, £~ and A, (c) the u-d-s quark matter with finite strange quark 
mass m s . 



easy to find the analytic solution of the above conditions; e.g. the proton frac- 
tion in a neutron star for a given baryon number density p reads p p /p n — 



1 + (ml/Sir 2 p n 



.2/3' 



-3/2 



This is a monotonically increasing function of p n 

and approaches to the asymptotic limit 1/8 from below. As the neutron density 
further increases and the electron chemical potential exceeds the muon mass, 
p c > m (1 = 105 MeV, the system composed of n, p, e~ and p~ is realized. 

As the baryon density increases further, hyperons enter into the game. This 
is because the Fermi energy of the neutron exceeds the threshold of the neutron- 
decay into hyperons. See Fig |25.6r b). Hyperons such as S _ and A may appear 
for p > (2 — 3)po- Which hyperon appears first depends on the still uncertain 
hyperon-nucleon interactions. 



25.3.2 Quark matter 

As the baryon number density p of the system exceeds (3 — 5)po, the neutrons, 
protons and hyperons start to percolate each other due to their finite sizes |27j . 
The quark number density p q for each flavor is related to p as p q = (N c /Nf)p 
with N c = 3 being the number of colors and Nf the active number of flavors so 
that the critical quark chemical potential, above which the percolation to quark 
matter takes place, is estimated as 

= (vr 2 /9g) V3 ^ (380 - 450) MeV. (25.3.8) 

Here we have assumed a non-interacting and degenerate quark matter composed 
of massless u and d quarks. 

Let us now consider the quark matter composed of only u, d and e _ . The 
condition of chemical equilibrium (d <-> u+e~), charge neutrality and the baryon 
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number conservation read ^d = Pu + Pc fPu — \p&~ Pc = 0, and |(p u +pd) = P- 
The factors 2/3 and —1/3 originate from the electric charges of the quarks and 
1/3 from the baryon number of a quark. If we assume non-interacting quarks at 
high density (fi q 3> m q ), one immediately finds fi u ~ 0.80 fid- Thus the Fermi 
energy of the d quark is slightly higher than that of the u quark, which is different 
from the situation of neutron matter where n and p have quite a different Fermi 
energies due to non-relativistic kinematics as shown in Fig l25.6f ay 

If the quark matter is composed of u, d, s and e _ , the chemical equilibration 
is achieved through the processes, d u + e~ , s^u + e - , and d + u <-> u + s. 
Then the equilibrium conditions read fid = fi u + /i c , fJ. s = /id, — \{Pd + ps) + 
|/5 U ~ p c = 0, and l(p u + Pd + Ps) = P- For fi q > m q , they lead to fi u = Pd = Ps 
and fj, e = 0. Namely the massless u-d-s quark matter is charge neutral by itself 
without electrons. If we have finite m s , then p s is reduced relative to p Uj d and 
the electrons become necessary to make the system charge neutral, as shown in 
FigHElc). 

25.4 Superfluidity in neutron-star matter 

If there exists attractive channel between the fermions near the Fermi surface, 
the system undergoes a transition to superfluidity or superconductivity in three 
spatial dimensions |28j . This is indeed the case in the neutron-star matter where 
the attraction between the protons due to spin-independent nuclear force in the 
(S, L, J) = (spin, orbital angular momentum, total angular momentum)=(0,0,0) 
channel leads to the condensation of 25+1 L,/= 1 So Cooper pairs (the proton 
superconductivity) and the attraction between the neutrons due to spin-orbit 
nuclear force in the (S, L, J)=(l,l,2) channel leads to the condensation of 3 P2 
Cooper pairs (the neutron superfluidity) [29 . 

Other than the superfluidity and superconductivity of the nucleons, the con- 
densation of pions (tt° and ir~) and that of kaons (K~) have been studied ex- 
tensively. For more details on nucleon superfluidity, meson condensation and 
its implication to the physics of compact stars, see the reviews [30l l3T| 132]. 

25.5 Color superconductivity in quark matter 

The quark matter exhibits the color superconductivity (CSC) which originates 
from the formation of Cooper pairs of quarks near the Fermi surface (see the 
recent reviews |33[ 134] and references therein) . Dominant attractive interaction 
responsible for the quark— quark pairing at high density is the color-magnetic 
interaction mediated by the gluon. 

There are some characteristic differences between CSC and the standard 
BCS-type superconductivity: 

(i) The quark matter at high density is a relativistic system where the quark 
chemical potential fi is comparable or larger than the quark mass m q . In 
such a case, the velocities of quarks near the Fermi surface are close to 
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light velocity and the magnetic interaction is not any more suppressed in 
comparison to the electric interaction. 

(ii) The color-magnetic interaction is screened only dynamically by Landau- 
damping, while the color-electric interaction is Debye screened as usual. 
Therefore, the collinear quark— quark scattering on the Fermi surface is 
dominated by the color-magnetic interaction, which leads to an uncon- 
ventional form of the fermion gap A k /ie" 

-c/^r s [35] Because of this 
non-BCS form where the coupling strength enters as ^foi^ instead of a s , 
A//x becomes sizable even in the weak coupling. 

(iii) Due to color and flavor indices of the quarks, the CSC gap acquires color- 
flavor matrix structure, which leads to various different phases depend- 
ing on n and T; they include the 2SC (2-flavor color superconducting) 
phase, the CFL (color-flavor locked) phase, the FFLO phase, the crys- 
talline phase, and so on as reviewed in [34] . 

25.5.1 The gap equation 

Let us illustrate the role of the above color-magnetic interaction by considering 
a simplified situation where quark matter is composed of only massless u and 
d quarks with equal Fermi energies at T — 0. In this case, the u quark with 
red-color and the d quark with green-color in flavor-singlet and color anti-triplet 
combinations are paired, while all the others quarks are unpaired [36j . This is 
called the 2SC phase. 

Using the standard Nambu-Gor'kov field = (g, <;*)', the Schwinger-Dyson 
equation for the quark self-energy X in the ladder approximation is written as 

m = ~ l j ^9 2 (p,k)rZS(p)r» b D a Jl(p-k), (25.5.9) 

where D"f v is the in-medium gluon propagator and is the bare quark-gluon 
vertex. The off-diagonal (anomalous) component of X in the Nambu-Gor'kov 
space is directly related to the gap function A(k) = A + (fc)A + + A_(/c)A_, 
where A± are the projection operators to the positive and negative energy 
states. Therefore, A + (k) and A_(fc) are interpreted as the quark gap and 
anti-quark gap, respectively. As for g{q,k), the Higashijima-Miransky ansatz 
(a momentum-dependent QCD coupling with phenomenological infrared regu- 
lator) may be adopted [37]. As for the gluon propagator, we take the screened 
propagator in the Landau gauge, 

/ ph P T \ 

DfJk) = 9 ^ - - -» ^ 8 a \ (25.5.10) 

' \ k 2 +m 2 D k 2 +ilml\k \/\k\) y ' 

where P^ T are the longitudinal and transverse projection operators. The lon- 
gitudinal (electric) part of the propagator has static screening by the Debye 
mass = (l/7r 2 )(7 2 /i 2 , so that the static interaction between quarks in the 
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Figure 25.7: (a) The quark gap A+(fe) as a function of |fe|//i for various quark 
chemical potentials, /j = 2" x A with A = 400 MeV and n = 1,2,3,12, ob- 
tained by solving the gap equation Eq. (|25.5.9[) numerically [35]. (b) Ratio of 
the coherence length £ c and the average inter-quark distance d q as a function of 
fi. 



coordinate space is the Yukawa-type short range potential. On the other hand, 
the transverse (magnetic) part has only dynamical screening due to Landau 
damping, so that the static interaction between quarks in the coordinate space 
is the Coulomb-type long range potential. 

In the high density limit (or equivalently the weak coupling limit due to 
asymptotic freedom) , only the interactions near the Fermi surface become rele- 
vant and the gap equation can be simplified to give the form [35j 

A+(|fc| = ft) S 2b n e -(3- 2 /V2)/ 3 (M) 5 (25.5.11) 



with b = 2567T 4 / g 5 (/z) and g(/z) = ^/Aira^n = fi) [34]. The characteristic form 
of the gap, e~ c / 9 , in Eq. (|25.5.1ip originates from the long-range color-magnetic 
interaction in Eq. (|25.5.10p and is different from the BCS form e~ c ^ 9 as we have 
mentioned. 

At low densities, sizable diffusion of the Fermi surface occurs, and the weak- 
coupling approximation in the asymptotic high density leading to Eq. (|25.5.1ip 
is not justified. Therefore, we solve the gap equation Eq. (|25.5.9p numerically 
with both magnetic and electric interactions [38]. In Fig l25.7f a). we show such 
numerical solution of A+(fc) as a function of |fe|//z for wide range of density. 
The figure shows that (i) the ratio of the gap and the quark chemical potential 
A + / /i can be sizable in magnitude of order 0.03 even for baryon density as high 
as p ~ lOpoj and that (ii) the gap function A+(fe) has non-trivial momentum 
dependence with a peak near the Fermi surface. 
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25.5.2 Tightly bound Cooper pairs 

The color superconductivity is a strongly coupled system partly due to the large 
value of coupling constant a s (n = fx) and partly due to long-range nature of the 
color-magnetic interaction. To clarify this point further, let us consider the 
coherence length £ c defined as a root-mean-square radius of the Cooper-pair 
wave function, <p + (r) oc (q(r)q(0)) oc f *fce«» ^ffij A . (fe) - 

Shown in Fig. 125. 7f b) is the ratio of £ c and the average inter-quark distance 
d q = (7r 2 /2) 1 / 3 ^~ 1 in the u-d quark matter. At high density, this ratio is very 
large (about 10 5 at fi ~ 10 6 MeV), while at low densities, the ratio becomes 
small (about 10 at /i = 800 MeV) and may even become less than 1 at lower 
densities. This situation is quite similar to the BCS-BEC crossover phenomenon 
[39] which was recently observed in ultracold atomic systems [40]: The result 
here suggests that the quark matter possibly realized in the core of neutron stars 
(fj, ~ 400 MeV) may be rather like the BEC of tightly bound Cooper pairs. For 
further studies of the BCS-BEC crossover in the relativistic system, see e.g. [4"T] . 



25.6 QCD phase structure 

The continuous QCD symmetry in the flavor-SU(3) chiral limit (m Uj d, s = 0), 
Eq. (l25.1.4|) . exhibits various symmetry breaking patterns depending on the tem- 
perature T and the quark chemical potential \x. In Table [25TT1 three examples of 
the symmetry realizations are shown; the QGP (quark-gluon plasma) phase, the 
NG (Nambu-Goldstone) phase and the CFL (color-flavor locked) phase. (CFL 
is one of the color superconducting phases to be discussed later in more details.) 
Each phase would appear in the T — /j, phase diagram as illustrated in Fig l25.8l 
In the intermediate values of /i ~ 400 MeV at low T, a variety of quantum 
phases have been proposed (30] [34] • 

As we have discussed in Sec l25. 1.121 the NG phase is characterized by the 
dynamical breaking of chiral symmetry due to nonvanishing chiral condensate 
(qq). On the other hand, in the color superconductivity at high chemical po- 
tential, the diquark condensate (qq) is formed as discussed in Sec l25.5l In the 
presence of these condensates, the light quarks and anti-quarks acquire the Dirac 
type mass M and the Majorana type mass A±, which leads to the relativistic 
quasi-particle spectrum of a quark near the Fermi surface, 

lo(p) = ^(^p 2 + M 2 - f i) 2 + \A+\ 2 . (25.6.12) 



25.6.1 Ginzburg-Landau potential for hot/dense QCD 

In this subsection, we focus our attention on the QCD phase structure at inter- 
mediate value of the chemical potential (/i ~ 400 MeV) where there would be 
an interplay between the quark— anti-quark pairing characterized by the chiral 
condensate (qq) and the quark— quark pairing characterized by the diquark con- 
densate (qq) [321 H3j ■ Study of this region is not only important to understand 
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Table 25.1: Symmetry breaking patterns of QCD in the flavor-SU(3) chiral limit 
(m u$ d,s = 0). The QGP phase, the NG phase, and the CFL phase denote the quark- 
gluon plasma phase, the Nambu-Goldstone phase, and the color-flavor locked phase, 
respectively. 



phase 


region 


unbroken continuous symmetries 


QGP phase 


"QCD 


SU(3) C <g> SU(3) L <g> SU(3) R ® U(1) B 


NG phase 


^- « 1, J^- « 1 

QCD QCD 


SU(3) ®SU(3) L+R ®U(1) B 


CFL phase 


3-£- < 1, jJi- » 1 

QCD QCD 


SU(3) C+L+R 



T 



~ 200 MeV 





QGP 




<qq> = 




<qq> - 


NG 




<qq> * 


y CFL 




T <qq> * 



~ 400 MeV 



Figure 25.8: The three basic phases of QCD in the T-fi plane for flavor-SU(3) chiral 
limit . 



the quantum phase transition to quark matter in the deep interior of the neu- 
tron stars, but also interesting in relation to similar phenomena in other systems 
such as the interplay between magnetically ordered phases and metallic super- 
conductivity [44j and that between superfluidity and magnetism in ultracold 
atoms [4"5] . 

To analyze such interplay in QCD in a model-independent manner, let us 
construct a Ginzburg-Landau (GL) potential f2 on the basis of the QCD symme- 
try Eq.(g5Xl as fX$A,d R ) = ^x($) + ^K + Here the 
chiral field <f>, which has a 3x3 matrix structure in the flavor space, is defined 
with the transformation property under Eq. (|25.1.3|) as, 

= (feJibJL), $ - e" 2 ** VL$Vl (25.6.13) 
On the other hand, the diquark field d L , which has a 3 x 3 matrix structure in 
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the flavor-color space, is defined with the transformation property as 



(25.6.14) 



where C = ij 2 ~f° is the charge conjugation matrix. Similar definition holds 
for e? R too. By definition, the 3x3 matrix [d L(R) ]ia belongs to the fundamental 
representation of SU(3) C and SU(3) L j R ). 

Most general form of the GL potential which is invariant under Q in Eq. (|25.1.4[) , 
written in terms of the chiral field up to C($ 4 ), reads [18], 



where "Tr" and "det" are taken over the flavor indices, i and j. The first three 
terms in the right hand side are invariant under Q <g> U(1) A , while the last term 
represents the axial anomaly which breaks U(1) A down to Z(6) A . The potential 
tt x is bounded from below for b\ + 62/3 > and 62 > 0. If these conditions are 
not satisfied, we need to introduce terms in C($ 6 ) to stabilize the potential, 
a situation we will indeed encounter. We assume cq to be positive so that 
the chiral condensate at low temperature is positive. Also, we assume that ao 
changes its sign at a certain temperature to drive the chiral phase transition. 

Most general form of the GL potential which is invariant under Q, written 
in terms of the diquark field up to 0(d 4 ), reads [4T>1 147j . 



n d = a Q TT[djl+d R d{} 

+(3, ([Tr(d L 4 )]2 + [Tr(d R <)] 2 ) + 2 (Tr[(djl ) 2 ] + Tr[(d R d R ) 2 ]) 



The transition from the normal state to color superconductivity is driven by ao 
changing sign. Unlike det <& in f2 x , terms such as det d L(R) are not allowed in 
Qd, since det rf L(R) carries baryon number and is not invariant under U(1)b- 

Finally, the interaction potential which is invariant under Q, written in terms 
of both chiral and diquark fields to fourth order, reads [4UI 138] , 



Q xd = 71 Tr[(d R 4)$ + (d L 4)$t] 

+Ai Tr[(d L 4)$$t + (d R dt )$t$] + x 2 Tr[d L 4 + d R d R ] ■ Tr[$t$] 



The term with the coefficient 71 originates from the axial anomaly which imposes 
the sign of 71 in Eq. (|25.6.17p and that of cq in Eq. (|25.6.15p being the same. 

Equations (|25.6.15[) . (|25.6.16p and (|25.6.17p constitute the most general form 
of the GL potential under the conditions that the phase transition is not strongly 
first order (i.e., the magnitudes of d L(R) are sufficiently smaller than those at 
zero temperature) and that the condensed phases are spatially homogeneous. 
To proceed analytically for the flavor-SU(3) chiral limit, we restrict ourselves to 
maximally symmetric condensates of the form: 




+/? 3 Tr[(d R 4)K4)]+/3 4 Tr(d L 4)Tr(d R 4). 



(25.6.16) 



+A 3 (det$ • TrKdjl)®- 1 } + h.c.) . 



(25.6.17) 



$ = diag(er, a, cr), d L 



d B = diag(c?, d, d), 



(25.6.18) 
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where a and d are assumed to be real and spatially uniform. We have chosen 
the relative sign between d L and d R in Eq. (|25.6.18p so that the ground state has 
positive parity, as is indeed favored by the axial anomaly together with finite 
quark masses. The above ansatz for the diquark condensate has residual sym- 
metry SU(3) C+L+R <g> Z(2) and is called the color-flavor locking (CFL) because 
of its symmetry realization |35]. (Note that Z(2) corresponds to the reflection, 

9l(r) > 9 L (R)-) 

The reduced GL potential with Eq. (|25.6.18p is 



f2 3 F = 



— < 

6 



V + ^ 
2 4 



~/d A a + \d A a 2 . (25.6.19) 



Here the axial anomaly leads to c > and 7 > 0, while microscopic calculation 
based on the Nambu— Jona-Lasinio model as well as the weak-coupling QCD 
suggests that A is positive and plays a minor role in comparison to 7 |48j . Note 
that we have introduced /-term (/ > 0) in case that b becomes negative. This 
system can have four phases with the following dynamical breaking patterns of 
continuous symmetries; 



QGP phase 
NG phase 
CFL phase 
COE phase 



a = 0,rf = 
a ^ 0,d = 
a = Q,d ^ 
a ^ 0,d ^ 



£^SU(3) C ®SU(3) L+R ®U(1) B 

g - su(3) c+L+R 

Q - SU(3) C+L+R . (25.6.20) 



The COE (coexistence) phase is favored by the axial anomaly, since the simul- 
taneous presence of d and positive a makes the GL potential lower because of 
the 7-term with 7 > 0. Note that even the unbroken discrete symmetry Z(2) is 
common between CFL and COE phases, so that they cannot be distinguishable 
from the symmetry point of view. 

For flavor-SU(2) chiral limit with m u ,d = and m s = 00, the condensates 
with the s quark disappear, so that we have $ = diag(er, a, 0) and rf L = — d R = 

diag(0,0, d). Then the reduced GL potential becomes Q2F = (f " 2 



4 U T 6 U 



\^d 2 + jd 4 J + \d l a l . Thus, the coexistence of d and a is disfavored in this 
case, because of the A-term with A > 0. 



25.6.2 Possible phase structure for realistic quark masses 

The mapping of the phase diagrams obtained from the GL potentials, O3F and 
in the a— a plane to the T—fi plane is a dynamical question which cannot be 
addressed within the phenomenological GL theory. Nevertheless, we can draw a 
speculative phase structure of QCD for m s ~ A QCD 3> m u .d 7^ by interpolating 
the phase structures obtained from f^p and f^F as shown in Fig l25.9l { 43] . In 
this figure, the double line indicates the first order phase transition driven by 
the negative b in Eq. (|25.6.19p . The single lines indicate the second order phase 
transitions (within the analysis of the GL potential without fluctuations) which 
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Figure 25.9: Schematic phase structure with two light (up and down) quarks 
and a medium heavy (strange) quark [43j . The double line indicates the first 
order transition. AY and HTYB are the second-order critical points at which 
the first-order line terminates. 



separate the d ^ and d = phases. We draw two critical points at which the 
first order phase transition turns into crossover; the one near the vertical axis 
indicated as "AY" (Asakawa-Yazaki critical point {50]) and the other one near 
the horizontal axis indicated as "HTYB" [43] . The latter is driven by the axial 
anomaly with positive 7 in Eq. (|25.6.19| . 

The existence of the AY critical point implies that the transition from the 
NG phase to the QGP phase on the \i — axis is a crossover. Indeed, the lattice 
QCD Monte Carlo simulations at finite T with the finite-size scaling analyses 
indicate that the thermal phase transition at fj, = is likely to be crossover [5T] . 
We note here that the AY critical point has special importance to the fluctuation 
observables in relativistic heavy-ion collisions [52] and the determination of its 
location is highly called for both theoretically and experimentally. 

On the other hand, the existence of the HTBY critical point implies that the 
haronic matter (characterized by a > d > and the quark matter characterized 
by d > a > are continuously connected with each other and both are classified 
into the COE phase. This is intimately related to the idea of hadron-quark 
continuity, i.e. smooth transition from the superfluid/superconducting hadronic 
matter to the superconducting quark matter 53 , 54J . Indeed, there are evidences 
of the continuity not only for the ground state but also for the excitation spectra: 
Typical example is the continuity of the flavor-octet vector mesons in hadronic 
matter at low n and the color-octet gluons in quark matter at high /x (53] [55] . 
Unfortunately, the lattice QCD simulations have difficulty to treat the matter 
with fi/T 3> 1 because of the severe sign problem originating from the complex 
fermion determinant in the presence of /z [56j . Therefore, the quantitative study 
in this region is still an open issue. 
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25.7 Simulating dense QCD with ultracold atoms 

Ultracold atomic systems and high density QCD matter, although differing by 
some twenty orders of magnitude in energy scales, share certain analogous phys- 
ical aspects, e.g., BEC-BCS crossovers [57]. Motivated by phenomenological 
studies of QCD that indicate a strong spin-singlet diquark correlation inside the 
nucleon [58], we focus here on modeling the transition from the 2-flavor quark 
matter at high density to the nuclear matter at low density in terms of a boson- 
fermion system, in which small size diquarks are the bosons, unpaired quarks the 
fcrmions, and the extended nucleons are regarded as composite boson-fermion 
particles |59] . This would be a starting point to understand the quantum phase 
transition at n ~ fi c between the hadronic superfluid discussed in Scc l25.4l and 
the color superconductivity discussed in Sec l25.5l 

Recent advances in atomic physics have made it possible indeed to realize 
boson-fermion mixture in the laboratory. In particular, tuning the atomic inter- 
action via a Feshbach resonance allows formation of heteronuclear molecules, as 
recently observed in a mixture of 87 Rb and K atomic vapors in a 3D optical 
lattice |60| . and in an optical dipole trap [61] . 

Let us start from a non-relativistic boson-fermion mixture with Hamiltonian 
density, 



where <j> is the boson and ip the fermion field. The two internal states of the 
fermions are labeled by spin indices a = {f , J.}. For simplicity, we consider an 
equally populated mixture of n bosons and n fermions with rif = = n/2. 

The bare boson-fermion coupling <7 b f is related to the renormalized coupling 
g h{ and to the s-wave scattering length a h{ by 



where £j(fc) = k 2 /2rrii (i — b, f ) is the single-particle kinetic energy, m R is the 
boson-fermion reduced mass, and A = n/(2ro) is a high momentum cutoff with 
tq being a typical atomic scale. We assume an attractive bare b-f interaction 
(<7bf < 0), tunable in magnitude, with A fixed, so that the scattering length 
a h{ can change sign: a bf — > <7bf?rt R /(27r) for small negative gu, while a bf — > 
ro for large negative gt,{- We keep the bare boson-boson and fermion- fermion 
interactions fixed and repulsive (<7bb > 0, <?g > 0) for the stability of this system. 

In the regime of the weak bare b-f coupling where the dimensionless param- 
eter rj = — l/(n 1 / 3 a bf ) is large and positive, the system at low temperature is 



n 



+ 





(25.7.21) 




(25.7.22) 
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Table 25.2: Correspondence between the boson-fermion mixture in ultracold atoms 
and the diquark-quark mixture in high density QCD. 



cold atoms 


dense QCD 


b (bosonic atom) 


d (diquark), 


fj,l (fermionic atom) 


q-t i (unpaired quark) 


N| ; j_ (boson-fermion molecule) 


(nucleon) 


b-f attraction 


gluonic attraction 


b-BEC 


2-flavor color superconductivity 


N-BCS 


nucleon superfluidity 



a weakly interacting mixture of BEC of the b-bosons (b-BEC) and degenerate 
f-fermions. The induced interaction through the density fluctuation of b-BEC 
may also lead to the pairing of f's (f-BCS). 

On the other hand, in the regime of strong bare b-f coupling where the r/ is 
large and negative, bound molecules of b-bosons and f-fermions called composite 
fermions, N = (bf), are formed with a kinetic mass m N = m b +m f . The s-wave 
scattering length of two N's of opposite spins can be estimated by the exchange 
of constituent b or f [59] . 

tTI 

" a bf . (25.7.23) 



2 ™ , 

This is the same in magnitude but is opposite in sign from the scattering length 
between difermion molecules due to different statistics. It can be shown that 
this result is the leading order term in an extension of the present model to 
large internal degrees of freedom. 

Eq. (|25.7.23|) implies that the low energy effective interaction between com- 
posite fermions in the spin-singlet channel is weakly attractive; the stronger the 
bare b-f attraction the weaker the N-N attraction. Such an effective interaction 
causes composite fermions to become BCS-paired (N-BCS) below a transition 
temperature, 

7 / o \ V'' 

Tc(N-BCS) = — f-J E N e* /(2kFa ™ ) ■ (25.7.24) 

where £n = fcp/2m N is the Fermi energy of the N. 

The above analyses for large |ry| suggest a possible phase structures of boson- 
fermion mixtures in the T — rj plane as shown in Fig. I25.10f a) . At intermediate 
bare b-f coupling (77 ~ 0) where a transition from the b-BEC phase to N-BCS 
takes place, the phase diagram would have complex structure depending on the 
relative magnitudes of gt>b, 9s, and g^. The f-BCS phase possibly occurs for 
77 > is not shown in this figure. For more detailed analyses of the phase 
diagram of the present model, see [55] • 
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Figure 25.10: (a) A possible phase structure of the boson-fermion mixture (such 
as 87 Rb and 40 K) in ultracold atoms with attractive b-f interaction and repulsive 
b-b and f-f interactions. Large and positive (large and negative) 77 corresponds 
to the weak (strong) b-f attraction, (b) A possible phase structure of QCD. 
Large (small) chemical potential fi corresponds to the weak (strong) coupling 
due to asymptotic freedom. 

The phase boundary in the region 77 ~ may be classified by the real- 
ization of internal symmetry. If we focus only on the continuous symmetries, 
the Hamiltonian density, Eq. (|25.7.21j) . has U(l) b ® U(l) f ® U(l) f symmetry 
corresponding to independent phase rotations of <f>, ipf and ip±. On the other 
hand, b-BEC and N-BCS break U(l) b and t/(l)b+(f T +f,) symmetries, respec- 
tively. The difference in such symmetry breaking patterns implies the existence 
of a well-defined phase boundary between b-BEC and N-BCS as indicated in 
Fig. 125.1'OT a). This is in contrast to the continuous BEC-BCS crossover in two- 
component Fermi systems. 

The phase structure we find for boson-fermion mixture of ultracold atoms 
displays features of that in QCD with equal numbers of u and d quarks. The 
ground state of such system at high density is the 2-flavor color superconduc- 
tivity (2SC) discussed in Sec l25.5.Tl The order parameter for color-symmetry 
breaking is the diquark condensate (d^) with the diquark operator d c = 
tij e abc [q] a C75 [q] \- The gap is of order a few tens of MeV; remaining quarks are 
unpaired and form degenerate Fermi seas. On the other hand, the ground state 
of QCD at low density is the nuclear matter with equal numbers of protons 
and neutrons denoted by A/"| j , a superfluid state with a pairing gap of a few 
MeV |32j ; the order parameter for the spontaneous breaking of baryon- number 
symmetry U(1) B is the six-quark condensate (Af^Af^) — ({d a [q^Y a )(d b [qiY b )) . If 
we model the nucleon, of radius r M ~ 0.86 fm, as a bound molecule of a diquark 
(of radius r d ~ 0.5 fm) and an unpaired quark, we can make the correspondence 
between boson-fermion mixture of cold atoms and the diquark-quark mixture 
in QCD as shown in Table HO 
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Such correspondence can be also found between the phase diagram of ultra- 
cold atoms in Fig |25.10T a) and that of dense QCD in Fig l25.107 b). In particular, 
the BCS-like superfluidity of composite fermion (N) with a small gap is a nat- 
ural consequence of the strong b-f attraction as shown in Eq. (|25.7.23|) . which 
may explain why the fermion gap in nucleon superfluidity is order of magni- 
tude smaller than the gap in BEC-like color superconductivity. It is thus quite 
interesting to carry out the experiments of boson-fermion mixture in ultracold 
atoms for wide range of the boson-fermion attraction. 

Note, however, that tuning the coupling strength at fixed density is not 
possible in dense QCD matter because of the running coupling a s (n = /i); 
furthermore, dynamical breaking of chiral symmetry and its interplay with the 
color superconductivity have an important role in the quantum phase transition 
in QCD as discussed in Sec l25.6.Tl With these reservations in mind, we sug- 
gest that fuller understanding, both theoretical and experimental, of the boson- 
fermion mixture [59] as well as a mixture of three species of atomic fermions 
[45l [62] can reveal properties of high density QCD not readily observable in 
laboratory experiments. 

25.8 Conclusion 

In this Chapter, we have discussed thermal and quantum phase transitions in 
QCD. The former is relevant to the physics of hot matter in early universe right 
after the big bang, while the latter is relevant to the physics of dense matter 
in the interiors of neutron stars and quark stars. There are three fundamental 
QCD phases in the T — /i plane: the NG phase, the QGP phase and the CSC 
phase (FigESU). 

We have shown an interesting possibility of hadron-quark continuity in which 
the superfluid hadronic phase and the CSC phase are continuously connected 
with each other at low temperature due to the QCD axial anomaly fFig |25.9f . 
We have also discussed that the existence of the nucleon superfluidity may be a 
logical consequence of the tightly bound diquarks interacting with the unpaired 
quarks in the CSC phase. Although such a system with fi/T ^ 1 is difficult 
to be treated in lattice QCD Monte Carlo simulations at present, a mixture of 
ultracold atoms with different masses, different statistics and different internal 
degrees of freedom would provide us with an exciting new tool to study the 
essential features of quantum phase transitions in dense QCD (Fig l25.10|) . 
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